Morphology Scaling of Drop Impact onto a Granular Layer 
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We investigate the impact of a free-falling water drop onto a granular layer. First, we constructed 
a phase diagram of crater shapes with two control parameters, impact speed and grain size. A low- 
speed impact makes a deeper cylindrical crater in a fluffy granular target. After high-speed impacts, 
we observed a convex bump higher than the initial surface level instead of a crater. The inner ring 
can be also observed in medium impact speed regime. Quantitatively, we found a scaling law for 
crater radius with a dimensionless number consisting of impact speed and density ratio between the 
bulk granular layer and water drop. This scaling demonstrates that the water drop deformation is 
crucial to understand the crater morphology. 



What happens when a drop of liquid impacts a granu- 
lar layer? Despite recent developments in the fundamen- 
tal physics of granular systems and fluid dynamics, this 
simple question remains unanswered. Much attention has 
been paid to dry granular materials [JJ. The mechanical 
properties of wet cohesive granular matter have been al- 
most the only aspect examined so far '4 13| . Recently, 
solid projectile impacts onto a granular layer have been 
well studied |4l4l4|- In contrast, impacts between a drop 
and a hard wall or a fluid pool have been examined ex- 
tensively |15l4l7|. Drop impact dynamics is related to 
many industrial applications, e.g., ink-jet printing, rapid 
spray cooling, and surface coating. Nevertheless, a drop 
impact to a granular layer has not been investigated until 
quite recently 18|. Here, we focus on the drop-granular 
impact. Our simple experiments serve as a starting point 
for exploring this phenomenon, which has great potential 
applicability to various fields such as planetary science, 
material science, civil engineering, and agriculture. For 
instance, the drop-granular impact may help to under- 
stand a geological scale impact in which a projectile is 
destroyed completely by the impact 19]. It may also re- 
late to the fossil rain drops which are small circular pit- 
like depressions in fine grain sediment and whose origin 
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is still a subject of controversy [2C 

We perform simple granular impact experiments with 
a free-falling water drop. When the drop impacts a gran- 
ular layer, it rebounds from the surface and then slowly 
sinks into the granular layer. Afterwards, a crater re- 
mains as evidence of the impact. Impact speed v is con- 
trolled by free-fall height h 2Mj which ranges in this 
experiment from 10 to 480 mm. The granular grains, 
which are commercial SiC abrasives, possess nonspheri- 
cal shapes and polydispersity and their grain size Dg is 
varied from 4, 8, 14, 20, or 50 /zm. A small vessel (30 
mm in diameter, 10 mm thick) is filled with grains by 
hand and used as a target. The five grain sizes can be 
grouped into three classes in terms of the packing frac- 



much smaller than that of random close packing. Thus, 
the granular layer includes numerous pores. Such a low 
packing fraction in 10'^-/im-size grains results from their 
irregular shape and polydispersity [23J, |2J| . The radius 
of a drop is fixed at R^ = 2.4 ± 0.2 mm. 

Figure [1] shows a typical sequence of a drop- 
granular impact taken by a high-speed camera (TAKEX 
FC350CL) at 210 fps. The first four images (Fig.^a-d)) 
display the impact moment. We observe the great defor- 
mation of the drop, which forms the crater rim. This rim 
forms within approximately 10 ms. After this initial im- 
pact stage, the drop rebounds and undergoes attenuated 
oscillation for a while. Then, the drop remains still on 
the surface as shown in Fig. [Tfe) . The oscillation settling 
time is 0.4 s. Finally, the drop penetrates very slowly 
into the granular layer and then forms a ring shape in- 
side the crater rim (Fig.[Tl^f)); this sinking process takes 
10 s. The time scales of early deformation and the later 
sinking are 10~^ s and 10^ s, respectively. Thus, the time 
scale expands by about three orders of magnitude. 

From an inspection of all video data, we find that the 
first expansive drop deformation time scale td does not 
depend on Dg and h. It is always td ~ 10~^ s. This time 
scale seems to be determined purely by the properties of 
the water drop. The surface tension of the water drop, 
7 = 7.2 X 10~^ N/m and its mass m^, = 5.8 x 10~^ kg 
yield the time scale t^ = \/mw /7 — 2.8 x 10~^ s. t-y 



should be the period of drop oscillation [25|, |26[. This 



estimate is consistent with the experimental result, i.e., 
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tion: 0.31 {Dg = 4, 
0.50 (Dg = 50 fi-ai). 
tion is about 10%. 



8 Aim), 0.44 {Dg = 14,20 fim), and 

The uncertainty of the packing frac- 

These packing fraction values are 
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The shape of the crater is not limited to that shown 
in Fig. [TJ It depends on the impact speed v and grain 
size Dg. We classify the crater shapes into four types, 
and draw a phase diagram of crater shapes (Fig. [5]). At 
low impact speed, clear splashing cannot be observed, 
and the drop sinks very gradually, leaving a cylindrical 
crater. It seems that the penetrating water drop com- 
presses the fiuffy granular layer mainly by capillary ef- 
fect. We call this a "sink crater". As the impact speed 
increases, the splashing of grains begins to emerge, and 
a "ring crater" is created (Fig. [1]). Between these two 
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FIG. 1: Typical sequence of a drop-granular impact {Dg = 4 fira and h — 160 mm). The transition from drop expansion to 
receding occurs between (c) and (d) . Drop deformation creates a crater rim at this early stage. The slow sinking of the drop 
forms an inner ring at a later stage (e,f). 
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FIG. 2: (color online). Phase diagram of crater shapes re- 
sulting from drop-granular impact. Twelve free-fall heights 
and five grain sizes are used in the experiments. In each ex- 
perimental condition, two independent runs are conducted; 
i.e., 120 impacts are examined in total. The presence of two 
symbols at a point indicates that two types are possible at 
that point. We have four types of crater shapes depending on 
Dg and h; sink, ring, flat, and bump. The characteristic h at 
which the qualitative crater shape changes is 100 mm. 



cases, a "flat crater" can be observed in the case of rel- 
atively small grains. In this type, the global structure 
of the crater is similar to that of the ring type, but the 
central region is rather flat. A "bump crater" appears at 
a high impact speed and large grain size. In this regime, 
the drop splits into smaller parts upon impact, and the 
largest part remains at the center. It absorbs surround- 
ing surface grains at the impact stage, and forms a convex 
bump by the sedimentation of collected grains during its 
penetration. Thus, drop impact can create even a convex 
shape as well as concave craters. 

To characterize the crater shape precisely, its surface 
structure is measured by a line laser displacement sen- 
sor system (KEYENCE LJG030) ^. Typical results 
are shown in Fig. [S] According to the phase diagram, 
h = 100 mm is a characteristic free-fall height that cor- 
responds roughly to the boundary of the sink, flat, and 
bump types. Therefore, we show the cases oi h = 80, 120, 
and 160 mm in Fig. El We also show ft, = 10 and 480 mm. 



wchich correspond to the slowest and the fastest impacts, 
respectively. 

From the surface depth map data (Fig. [31Ja-c)), we 
compute the radial depth function d{r), where r is the 
distance from the center of the crater, and plot the cor- 
responding d{r) curves (Fig. [31[d-f)). When Dg = A ^m, 
the shortest height {h — 10 mm) impact creates a 
deeper crater than medium-height impacts oi h = 80 and 
120 mm. In medium-height impacts, the impact inertia is 
inadequate to make a deep crater, but it compresses the 
granular layer a little. And it effectively suppresses the 
capillary based compression. Therefore a lower impact 
energy is better able to deform a fluffy granular layer. 
The crater depth is not a monotonic function; when 
the impact speed is large enough, craters become deeper 
again (Fig. [3][a,d)). With large grains, the crater shapes 
are completely different as shown in Fig.[3l[b,c,e,f). They 
are no longer simple concave shapes as mentioned be- 
fore. Interestingly, the bump height is much higher than 
the initial granular layer level z — (Fig. Eff), for 
h = 480 mm). 

A solid projectile impact cratering was also measured 
using a laser profilometry [ll| . It was found that the 
crater shapes are very close to hyperbolic function. How- 
ever, the drop-granular impact creates very different 
craters due to the drop deformation and capillary effect, 
as mentioned above. 

Next, we analyze the characteristic length scales. First, 
we discuss the vertical length scale. The minimum crater 
depth dmin with respect to free-fall height h is shown in 
Fig. mja). The dmin curves for small Dg have a peak at 
a certain h (^ 100 mm). On the other hand, the dmin 
curves for large Dg are rather weakly decreasing functions 
oih. 

Analysis of the horizontal length scale R, the radius 
of the crater, is more straightforward. Here, we assume 
R is determined by the deformed drop radius Rd- This 
assumption is roughly confirmed in Fig. [T] (c,d). Water 
drop deformation by the impact is estimated by Oku- 
mura et al. [25|, |26|. They studied water drop impact 
onto a super- hydrophobic substrate [28|. Using Euler's 
equation, the Laplace pressure gradient, and geometrical 
conditions, they have derived a horizontal deformation 
Rd that is scaled as Rd ~ R^We^^'^, where We is the 
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FIG. 3: (color online). Grayscale surface depth maps of Dg = (a) 4, (b) 14, and (c) 50 ^m cases. The vertical height z = \s 
the initial granular surface level before impact. All depth map data are shown in 15 x 15 mm'^ squares with axes in units of 
mm. The corresponding radial depth d as a function of distance r from the center of the crater is shown in (d-f). As (a,d) show, 
ft = 10 mm results in deeper craters than /i = 80 or 120 mm cases. In large Dg impacts, it is hard to observe a hemispherical 
crater. Furthermore, the central bump is higher than the initial height 2 = when Dg = 50 /xm and h = 480 mm. This is 
called a bump crater in the phase diagram shown in Fig. (2] 



Weber number. We — 2p^i?u,w^/7, and pw is the density 
of water. We is proportional to h through the free-fall 
relation w^ ~ h. Thus, R should scale as -R '^ h^/^. 
The actual R data and this scaling slope are shown in 
Fig. Hljb). Although the global structure seems to be 
sound, the data spread depending on Dg. We believe this 
spreading is caused by the difference between the density 
of the bulk granular layer pg and that of the water drop 
Pui- Of course, pg is essentially the same parameter as 
the packing fraction. If the ratio Pgj Pw is large, the drop 
is greatly deformed by the impact. Consequently, a large 
crater is created. The pg is the simplest representative 
property of bulk granular systems. Thus, it is natural 
to consider the bulk density ratio when discussing the 
granular impact physics. In fact, this bulk density ratio 
is also useful for discussing the dynamic scaling of solid 
projectile impacts to a granular layer [29|. Using this 
density ratio, we finally define the scaling as. 
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To check this scaling, the normalized crater radius 
R/Ryj is plotted as a function of the scaling param- 
eter [pgl pyj)We^/'^ (Fig. m^c)). This scaling implies 
R ~ Pg{R^v'^ I ptolY^^ ■ As expected, the data collapse 
to a line. Although the original scaling Rd ~ We^/^ was 



deduced for a hard wall target, it is still valid for a de- 
formable granular target. At large We, data deviate from 
the scaling due to the drop splitting by the high impact 
energy. 

Similar crater radius scaling was found in solid projec- 
tile impact experiments [J, la, [7[ . The scaling exponent 
1/4 means that the most of the impact energy is spent 
in lifting the ejecta resulting from cratering [J|. If this 
is true, the crater size should be smaller with larger pg. 
However, the current result shows the opposite tendency. 
The physics of solid impacts and drop impacts are quite 
different, even though the crater radius scaling exponent 
is similar. In drop-granular impact, drop deformation 
effect is much more crucial than ejecta splashing. The 
impact Reynolds number / = We^/^Re^^"^ is derived by 
substituting the water drop characteristic length scale 
and time scale into the definition of the Reynolds num- 



ber Re 



/vt-y/v with kinematic viscosity v []j 



systematic evaluation of the impacting fluid's surface ten- 
sion and viscosity dependency is necessary for using an 
advanced dimensionless number such as /. This is an 
open problem. 

We have demonstrated a simple but original experi- 
ment of drop-granular impact. The drop deformation 
time scale can be evaluated by the surface-tension-based 
elasticity of the water drop. We classified the result- 
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FIG. 4: (color online). Characteristic length scales, (a) Mini- 
mum crater depth dmin and (b) crater radius R, as functions 
of free-fall height h. dmin shows a peak at /i ~ 100 mm at 
Dg = 4 and 8 /im cases (a). Other curves are slightly de- 
creasing functions of h. R roughly agrees with the scaling 
R ~ h^' (b). To collapse all the R data, a dimensionless plot 
R/Rui vs {pg/ pm)We ' is shown (c). We can confirm good 
data collapse to the scaling. In all plots, symbols and colors 
indicate grain size Dg of the granular layer, as shown in the 
legend. 



ing crater shapes into four characteristic types and com- 
pleted the phase diagram. We found that the crater 
radius scales as R/R^ ^ [Pg/ Pw)We^^^. However, the 
vertical length behaves in a more complex manner. It 
seems to be affected by competition between multiple 
fluid-granular interactions. To reveal the effect more ac- 
curately and obtain more universal scaling, we must sys- 
tematically vary many parameters, e.g., viscosity, sur- 
face tension, size of impacting drop, and grain shape and 
surface properties. This report is the first step toward 
revealing this rich drop-granular impact world. 
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